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Abstract 

We provide a scenario in which a hidden £7(1) gauge boson constitutes dark matter of the 
Universe and decays into the standard-model particles through a kinetic mixing with an U(1)b-l 
gauge boson. Interestingly, our model can naturally account for the steep rise in the positron 
fraction recently reported by PAMELA. Moreover, we find that due to the charge assignment of 
U{V)b-Li only a small amount of antiprotons are produced in the decay, which is also consistent 
with the PAMELA and other observational data. 
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I. INTRODUCTION 



The presence of the dark matter has been firmly established by numerous observations, 
although we have not yet understood what dark matter is made of. It is promising that the 



PAMELA [lj] and Fermi (formerly GLAST) '. 
of dark matter, helping us to pin down the dark matter particle. 

Recently much attention has been attracted to a scenario that the dark matter decays 
into the sta.dard-.ode! (SM, pa rtlcleS g fl fl fl fl Q fl fl Q, S1 oee the eoergetie par- 
tides produced in the decay may account for the excesses in the positron fraction reported 



by PAMELA 12J or HEAT 13| and in the diffuse extra galactic gamma-ray background 
observed by EGRET Q. The PAMELA also released data on the antiproton flux 16], 
suggesting that most of the observed antiprotons are the secondaries produced by interac- 
tions between the primary cosmic rays and the interstellar medium. The suppression in the 
antiproton flux from the dark matter decay is particularly important because some decaying 
dark matter scenarios predict too large antiproton flux at the solar system 6]. 

We have recently proposed a scenario that a hidden U(l) gauge boson may decay into 
the SM particles through its kinetic mixing with the U(l)y in the SM gauge group [8j. In 
order to explain the longevity of the hidden gauge boson, the kinetic mixing needs to be 
suppressed down to an extremely small value, and direct couplings between the hidden and 
SM sectors must be negligibly small 17|. In order to achieve those features we introduced 
messenger fields having a large mass close to the grand unification theory (GUT) scale, and 
assumed that the hidden sector is sequestered from the SM sector. 

In this letter we present a model in which such a small coupling between the hidden U(l) 
gauge boson and the SM particles are naturally realized in the extra dimension framework. 
The essential idea is to separate the hidden sector from the SM sector in an extra dimension 
so that the two sectors can communicate with each other only through the interactions with 
another U(l) m gauge field in the bulk, which is assumed to be broken at a high-energy 
scale. The branching ratios simply reflect the charges of the SM particles under the U(l) m 
gauge group. Interestingly, our model can account for the steep rise in the positron fraction 
reported by PAMELA as well as the gamma-ray excess observed by EGRET. Moreover, the 
antiproton flux can be suppressed enough to be consistent with the PAMELA observations, 
if U(l) m is identified with U(1)b-l- 
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II. SET-UP 



We introduce an extra dimension with two branes at the boundaries. Suppose that the 
hidden gauge sector is on one brane and the SM particles are on the other brane well 
separated from each other. In such a set-up, dangerous direct interactions between the two 
sectors are exponentially suppressed, and the hidden gauge boson becomes quasi-stable * 1 . 
It is worth noting that its stability is guaranteed by the geometric separation in the extra 
dimension, not by a discrete symmetry such as a .R-parity. If there is a U(l) b-l gauge field in 
the bulk, the hidden U(1)h gauge field can have an unsuppressed gauge kinetic mixing with 
the U(1)b-l- After integrating out the heavy U(1)b-l gauge boson, the effective couplings 
between the hidden U(1)h gauge boson A H and the SM particles are induced, which enables 
Ah to decay into the SM particles. The longevity of the hidden gauge boson is realized by 
the hierarchy between B — L symmetry breaking scale and the weak scale. As we will see 
below, taking the B — L breaking scale around the GUT scale * 2 and the mass of the hidden 
gauge boson of (9(100) GeV naturally leads to the lifetime of (9(10 26 ) second that is needed 
to account for the positron excess. 

Let us first consider the kinetic mixing between the hidden U(1)h and the U(1)b-l gauge 
symmetries. The relevant effective interactions in the four dimensions are written as 

+ l -ntA Hll A H + l -M 2 A B ,A% (1) 

where A denotes a coefficient of the kinetic mixing of order unity, and F^ and F^ are the 
field strengths of the U(1)h and U (1)b-l gauge bosons, Ah and Ab, respectively. We assume 
that both gauge symmetries are spontaneously broken, and therefore Ah and Ab acquire 
masses m and M, respectively. We also assume that the kinetic mixing is unsuppressed as 



The interactions are suppressed by e~ M * L , where M* is the five-dimensional Planck scale and L denotes 
the size of the extra dimension. For M*X ~ 10 2 , the direct couplings are so small that the hidden 
gauge boson will become practically stable in a cosmological time scale. Then is roughly equal to 
Alp/ 10 ~ 10 17 GeV which is larger than the GUT scale (~ 10 15 GeV), and so, our analysis in the text is 
valid. 

# 2 The seesaw mechanism 18j for neutrino mass generation suggests the mass of the right-handed neutrinos 
at the GUT scale ~ 10 15 GeV. Recall that the right-handed neutrinos acquire the masses from the U(1)b-l 
gauge symmetry breaking. Thus, it is quite natural to consider the B — L breaking scale around the GUT 
scale of ~ 10 15 GcV. 



3 





quark 


lepton 


N^B-L) 2 


l 
3 


1 



TABLE I: The coefficients appearing in the partial decay rates for quarks and leptons. 

A = 0(0.1), and we take m = 0(100) GeV and M = 0(1O 15 ) GeV throughout this letter. 
We can make the kinetic terms canonical and diagonalize the mass matrix by appropriate 
transformations. The relations between (Ah, Ab) and the mass eigenstates (A' H ,A' B ) are 

2 \ 
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where we have approximated m 2 <C M 2 and A ^5 0.1 for simplicity. 

The low-energy effective interactions between the hidden gauge boson A' H and the SM 
fermion ipi can be extracted from the U(1)b-l gauge interactions by using the relation Q, 



772 

£ int = qiA^ipa^ipi D -Xqi—A^ipij^i 



(4) 



where denotes the B — L charge of the fermion ip^. The partial decay width for the SM 
fermion pair is 



H 



m \ 4 
— m, 



(5) 



12tt \MJ 

where we have neglected the fermion mass, and N{ is the color factor (3 for quarks and 1 
for leptons). Thus the lifetime r is given by 

2.5 x 10 27 sec / m \ ^ ( M 



T 



(6) 



V^Nitf VIOOGeV/ V!0 15 GeV 
where the sum is taken over those SM fermions of masses lighter than m/2. 

We show the coefficient A^g 2 for the quarks and the leptons in Table [B It should be noted 
that the branching ratios are not sensitive to the mass of Ah and they simply reflect the 
B — L charge assignment, which makes our analysis very predictive. 



III. COSMIC-RAY SPECTRA 



In this section we show the predicted spectra for the positron fraction, gamma-ray and 
antiproton fluxes based on the decay modes shown in the previous section. More precisely, 
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FIG. 1: The predicted positron fraction from Ah decay via the kinetic mixing with U(1)b— l (blue 
line) and f (1)5 (magenta line) using the M2 propagation model, compared with the experimental 
data, including the recent PAMELA results. 

the branching ratios are 2/39(2/37), 2/13(6/37) and 1/13(3/37) for a quark pair, a charged 
lepton pair and a light neutrino pair, respectively, if the top quark decay channel is (not) 
allowed kinematically. To estimate the spectra of gamma, positron and antiproton, we use 



the PYTHIA 



191 ] Monte Carlo program. After cosmic-ray particles are produced during the 



decay of Ah, the following calculations are straightforward and identical to those adopted in 



Rcf. 



and so, we show only the final results in this letter. 



or readers who are interested 



in the details of the calculations should be referred to Ref. |6| and references therein. 

In our numerical calculations we set m = 300 GeV and the lifetime r = 1.2 x 10 26 
seconds, and we use the so-called M2 diffusion model that are consistent with the Boron 
to Carbon ratio (B/C) and produce the m in lma > positron nnx fl Q. In Fi, ffl (Hue 
line), we show the predicted positron fraction together with the recent PAMELA data and 
other experiments. The prediction of our model fits very well with the excess reported by 
PAMELA. The positron fraction steeply increases from E ~ 10 GeV and drops off sharply 
at E = m/2 = 150GeV, which is mainly caused by the contribution of e + directly produced 
by the Ah decay. Such a drop-off can be checked by the upco ming PAMELA data in the 
higher energy region. For the MED and Ml diffusion models [20j, the positron fraction 
becomes slightly softer in the low energy (around several tens GeV), while the difference 
is negligible in the high end around the peak. We mention here that, with these new data 



of positron from PAMELA 



12j, the background estimation may be different from what we 



adopted here 2l|, [22|. However, since the signal of positron from decay of Ah is negligible 
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FIG. 2: The predicted gamma-ray flux of Ah decay via the kinetic mixing with U(1)b-l (blue 
line) and U(l)s (magenta line), compared with the EGRET data 



in the low energy and is important in the high energy region, i.e. E > 10 GeV, we expect 
that our model will still be able to explain the excess even with such a new background 
estimation that fits the low energy data better. 

The gamma-rays are mainly produced by the 7r° generated in the QCD hadronization 
process, since quark pairs are produced from the decay of Ah- In Fig. [2] (blue line), we plot 
the gamma-rays together with the EGRET data. The gamma-ray flux peak at E ~ 20 GeV, 
and ends at E = m/2, which reflects the mass of decaying dark matter. With no surprise, 
we see that the excess observed by EGRET may also be explained by the decay of Ah ■ 

Finally we show in Fig. [3] our predicted contribution to the antiproton flux (blue line). In 
Fig. [3l we have not included the prediction on the secondary antiproton flux, which should 
explain the BESS data 23j. Importantly, the predicted contribution to the antiproton flux 
from the Ah decay is smaller than the observed one by more than one order of magnitude, if 
the MIN propagation model 24J is adopted, as shown in Fig. [3] (a). The suppression in the 
antiproton flux is particularly crucial, because the recent data from PAMELA are consistent 
with the previous BESS result, which suggested that the secondary production dominates 
the observed antiproton flux. Furthermore, too many antiprotons tend to be generated as a 
by-product when we require the dark matter annihilation/decay to account for the positron 
excess f], Of course, the predicted antiproton flux still has a large uncertainty mainly 
due to our poor understanding of the cosmic-ray propagation inside our galaxy. As we can 
see in Fig.|3](b), for different propagation models (MED and MAX), the antiproton flux from 
decaying Ah can be enhanced by about two orders of magnitude. However, our scenario 
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FIG. 3: (a) Antiproton flux from the Afj decay via the kinetic mixing with U(1)b—l (blue line) 
and U(l)5 (magenta line), compared with the BESS data, (b) Predicted antiproton fluxes from the 
Ah in the U(1)b—l case, using different propagation models. Note that the secondary antiproton 
contributions are not included. 

can still be consistent with the observed antiproton flux.* 3 
IV. DISCUSSION AND CONCLUSIONS 



So far we have focused on the case that the U(l) m gauge symmetry in the bulk is identified 
with the U(1)b-l symmetry, but it is also possible to consider another anomaly-free charge 
assignments given by a linear combination of the B — L and the hypercharge Y. However, as 
we can check easily, the hadronic decay branching ratio, which is a measure for the antiproton 
flux, becomes smaller as the contribution of the hypercharge is suppressed. In this respect, 
our choice of the U(1)b-l is well motivated by observation, since the antiproton flux is 
naturally suppressed enough to be consistent with the PAMELA data. For comparison, let 
us consider a Z7 (l)s, so-called "fiveness", instead of U(1)b-l- The charge Q$ of the U (1)5 is 
given by the following linear combination of the charges under U(1)b-l an d U{l)y 



2V10 



Y-(B-L). 



(7) 



# 3 We notice that, the estimated background secondary antiproton flux drops quickly as the energy increases, 
and that the signal for the MED diffusion model seems to exceed the PAMELA data on the antiproton 
to proton ratio in the high energy region. However, our scenario is consistent with the PAMELA data in 
the high end if the MIN propagation model is adopted. 
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We have similarly estimated the spectra for the positron fraction, gamma-rays, and antipro- 
tons when the U(l) m in the bulk is identified with {7(1)5. In order to compare with the 
U(1)b-l case, we take the same parameter sets, i.e. m = 300 GeV and r = 1.2 x 10 26 sec. 
The decay branching ratios of Ah decay are roughly 4.5%, 18.3%, 6.5%, 4.8% and 2.3% 
for up-, down-type quark, charged lepton, neutrino and W-boson pairs, respectively. We 
show the results of cosmic rays as the magenta lines in Figs. [U E] and [31 The distribution of 
positron fraction in fiveness case seems flatter compared to the B-L case and the inclination 
of the turnup is softer. This is due to the fact that positrons generated from the decay of 
hadrons (mainly from ir + ) become more important and the branching ratio of e + e~ mode 
decreases simultaneously. The fiveness model predicts a larger excess in the gamma-ray flux, 
since more ir°s are produced. Finally, we can see from Fig. [3] that the antiproton flux is 
enhanced for the C/(l)g compared to the case of U(1)b-l due to a larger total branching 
ratio of quark pairs. Note that, in order to explain the PAMELA positron fraction data in 
the high energy region, we need to adopt a shorter lifetime for fiveness model compared to 
the case of B — L. Accordingly, the gamma-ray and antiproton fluxes increase, which may 
result in a tension with observed data and make the fiveness model being slightly disfavored. 



In our previous work jg] we considered a scenario that the hidden U(1)h gauge boson mixes 
with the {7(1) y- The predicted spectra in this case are somewhat between the B — L and 
fiveness cases. 

Let us comment on the production of the hidden gauge boson in the early Universe. 
Although the couplings of Ah to the SM particles are extremely suppressed, we can generate 
a right abundance of Ah from thermal scatterings as follows. For the reheating temperature 
about 10 15 GeV, the B — L gauge bosons will be in thermal equilibrium, and the Ah can be 
produced through a non-renormalizable coupling such as 

where M* denotes the five-dimensional Planck scale (see the footnote (1)) and k is a nu- 
merical coefficient of order unity. The presence of such non-renormalizable operator on the 
hidden brane is natural since it is allowed by the gauge symmetries. The abundance of Ah 
produced via the operator is roughly estimated by 

n „*» = 0( o.i) (^) 2 (jy 3 (9 ) 



where Tr is the reheating temperature. Note that the required reheating temperature is 
smaller than the cut-off scale M*. Also non-thermal production of Ah by e.g. the inflaton 



decay 
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291 ] should work as well. 
In this letter, we propose a scenario that a hidden gauge boson constitutes the dark 
matter of the Universe and decays into the SM particles via the kinetic mixing with the 
U(1)b-l gauge field in the bulk. Interestingly, our model can account for the steep rise in 
the positron fraction reported by PAMELA as well as the gamma-ray excess seen by EGRET, 
while avoiding the constraint on the antiproton flux by PAMELA and other experiments, 
if the M2 propagation model is adopted, due to the smallness of quark's quantum number 
under the gauged U(1)b-l- Moreover, the very small decay rate of the hidden Ah gauge 
boson dark matter is realized naturally by the hierarchy between weak scale and the large 
B — L breaking scale which is about the GUT scale as suggested by the neutrino masses. 

Acknowledgments 

This work was supported by World Premier International Research Center Initiative 
(WPI Initiative), MEXT, Japan. 



[1] P. Picozza et al, Astropart. Phys. 27, 296 (2007) [arXiv:astro-ph/0608697~] ; see also the 



webpage: http:/ /pamela.roma2. infn.it/index.php 



[2] Fermi Gamma-ray Space Telescope (formerly GLAST) collaboration, see the webpage: 



http:/ /fermi. gsfc.nasa.gov/ 



[3] F. Takayama and M. Yamaguchi, Phys. Lett. B 485, 388 (2000) [arXiv:hep-ph/0005214] . 
[4] W. Buchmuller, L. Covi, K. Hamaguchi, A. Ibarra and T. Yanagida, JHEP 0703, 037 (2007) 
[arXiv:hep-ph/0702184] . 

[5] G. Bertone, W. Buchmuller, L. Covi and A. Ibarra, JCAP 0711, 003 (2007) [arXiv: 0709 .22991 
[astro-ph]]. 

[6] A. Ibarra and D. Tran, Phys. Rev. Lett. 100, 061301 (2008) [arXiv: 0709 .45931 [astro-ph]]; 

JCAP 0807, 002 (2008) [arXiv:0804.4596l [astro-ph]]. 
[7] K. Ishiwata, S. Matsumoto and T. Moroi, arXiv:0805.lT33l [hep-ph]. 



9 



[8] C. R. Chen, F. Takahashi and T. T. Yanagida, larXiv:0809.0792l [hep-ph]. 
[9] A. Ibarra, A. Ringwald and C. Weniger, larXiv:0809.3l96l [hep-ph]. 
[10] C. R. Chen and F. Takahashi. larXiv:0810.4lT0l [hep-ph]. 

[11] P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu. larXiv:0811.0176l [hep-ph], 
[12] O. Adriani et al. larXiv:0810.4995l [astro-ph], 

[13] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) 

|arXiv:astro-ph/9703l"92] . 
[14] P. Sreekumar et al. [EGRET Collaboration], Astrophys. J. 494, 523 (1998) 

[arXiv:astro-ph/9709~257] . 
[15] A. W. Strong, I. V. Moskalenko and O. Reimer, Astrophys. J. 613, 956 (2004) 

|arXiv:astro-ph/0405441| . 
[16] O. Adriani et al. larXiv:0810.4994l [astro-ph]. 

[17] M. Pospelov, A. Ritz and M. B. Voloshin, Phys. Lett. B 662, 53 (2008) [arXiv:0711.4866l 
[hep-ph]]; J. L. Feng and J. Kumar. larXiv:0 803.4196 [hep-ph]; D. Hooper and K. M. Zurek, 
Phys. Rev. D 77, 087302 (2008) |arXiv: 080 1.36861 [hep-ph]]; J. March-Russell, S. M. West, 
D. Cumberbatch and D. Hooper, JHEP 0807, 058 (2008) |arXiv:0801.3440l [hep-ph]]. 

[18] T. Yangida, in Proceedings of the "Workshop on the Unified Theory and the Baryon Number 
in the Universe", Tsukuba, Japan, Feb. 13-14, 1979, edited by O. Sawada and A. Sugamoto, 
KEK report KEK-79-18, p. 95, and "Horizontal Symmetry And Masses Of Neutrinos" , Prog. 
Theor. Phys. 64 (1980) 1103; M. Gell-Mann, P. Ramond and R. Slansky, in "Supergravity" 
(North-Holland, Amsterdam, 1979) eds. D. Z. Freedom and P. van Nieuwenhuizen, Print-80- 
0576 (CERN); see also P. Minkowski, Phys. Lett. B 67, 421 (1977). 

[19] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) [arXiv:hep-ph/0603175| . 

[20] T. Delahaye, R. Lineros, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 77, 063527 
(2008) |arXiv:0712.23T2l [astro-ph]]. 

[21] I. V. Moskalenko and A. W. Strong, Astrophys. J. 493, 694 (1998) [arXiv:astro-ph/9710124~] . 

[22] E. A. Baltz and J. Edsjo, Phys. Rev. D 59, 023511 (1999) |arXiv:astro-ph/9808243| . 

[23] S. Orito et al. [BESS Collaboration], Phys. Rev. Lett. 84, 1078 (2000) 
[arXiv:astro-ph/9906426] . 

[24] F. Donato, N. Fornengo, D. Maurin and P. Salati, Phys. Rev. D 69, 063501 (2004) 
|arXiv:astro-ph/0306207] . 

10 



[25] M. Cirelli, M. Kadastik, M. Raidal and A. Strumia. larXiv:0809.2409l [hep-ph]. 
[26] M. Fujii and K. Suzuki, Phys. Rev. D 66, 035002 (2002) [arXiv:hep-ph/0202038l . 
[27] M. Endo, M. Kawasaki, F. Takahashi and T. T. Yanagida, Phys. Lett. B 642, 518 (2006) 
[arXiv:hep-ph/0607170] . 

[28] M. Endo, K. Hamaguchi and F. Takahashi, Phys. Rev. D 74, 023531 (2006) 
|arXiv:hep-ph/06050"9T] . 

[29] M. Endo, F. Takahashi and T. T. Yanagida, Phys. Lett. B 658, 236 (2008) 
[arXiv:hep-ph/0701042] ; Phys. Rev. D 76, 083509 (2007) [ arXiv: 0706 .09861 [hep-ph]]. 



11 



